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ABSTRACT 

The NASA K u -band Radar Ocean Wave Spectrometer (ROWS) is an experimental proto- 
type of a possible future satellite instrument for low data rate global waves 
measurements. The ROWS technique, which utilizes short-pulse radar altimeters 
in a conical scan mode near vertical incidence to map the directional slope 
spectrum in wave number and azimuth, is briefly described. The potential of the 
technique is illustrated by some specific case studies of wave physical processes 
utilizing the aircraft ROWS data. These include i) an evaluation of numerical 
hlndcast model performance in storm sea conditions, 1i) a study of fetch-limited 
wave growth, and 1 i i ) a study of the fully-developed sea state. Results of 
these studies, which are briefly summarized, show how directional wave spectral 
observations from a mobile platform can contribute enormously to our understand- 
ing of wave physical processes. 

1. INTRODUCTION 

A simple method for measurmq the vector wave number spectrum of ocean sur- 
face gravity waves from aircraft and satellite platforms using modified radar 
altimeters has been described and investigated theoretically by Jackson (1981) 
and demonstrated experimentally by Jackson et aK (1984a). In this paper, rather 
than dwell on the details of the technique, we will present some specific air- 
craft results that will serve to Illustrate the enormous potential of this meas- 
urement technique for furthering our understanding of wave physical processes. 

The GSFC K u -band Radar Ocean Wave Spectrometer (ROWS) Is a noncoherent, 
short pulse radar that uses a near-nadir directed conically scanning antenna to 
map wave directionality. Table 1 gives the pertinent Instrument characteristics. 
Figure 1 depicts the aircraft measurement geometry. A small rotarty antenna bore- 
sighted to 16° Incidence produces a footprint at the nominal 10 km aircraft al- 
titude measuring approximately 1500 m in the range (x) dimension and 700 m In the 
orthogonal a’Huthal (y) dimension. The surface Is probed In the range dimension 
using 12.5 ns compressed pulses. At the nominal Incidence angle for peak power 
return, e » 13°, the surface, range resolution Is 8 m. The directional resolution 
of the ROWS Is obtained by a simple phase front matching condition between elec- 
tromagnetic and ocean wave Fourier contrast wave components across the relatively 
broad azimuth beamwldth. That is, the broad beamwldth functions to Isolate, or 
resolve, ocean Fourier components whose wave vectors k^ * (k, ♦) are aligned with 
with the radar azimuth + . In this respect, the ROWS technique Is similar to the 
dual frequency technique Investigated by Alpers and Hasselmann (1978). Let the 
fractional cross section variation for any pixel (x, y) be denoted 6o/n. The 
fractional reflectivity modulation m seen by the radar Is given by Wcr aver- 
aged laterally across the beam. If G(y) denotes the lateral gain pattern, then 
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RADAR OCEAN WAVE SPECTROMETER 




Figure 1. ROWS Measurement Geometry, 


TABLE 1. ROWS INSTRUMENT CHARACTERISTICS 


Frequency: 

13.9 TTHz 

Pulse type: 

Linear FM. 100 MHz, 1.2 us chirp 

Pulse length: 

12.5 ns compressed 

Peak power: 

2 kW 

PRF: 

100 Hz max 

Detection : 

Noncoherent, square law 

Antenna : 

10° elevation X 4 3 azimuth printed 
circuit, 16° incidence boresiqht, 

6 rpm rotation rate 

Data : 

Digital, max 1024 six bit word frame 
size, sample gate width selectable 
2, 5, ... ns: recording at full PRF 


m(x, <f> )* <5W/W * f G 2 ^ |6g/q)dy (1) 

where 5W/W is the fractional modulation in received power (averaged over the 
clutter fluctuations). 

In the near-vertical specular backscatter regime In which the ROWS operates, 
the cross-section variation Is primarily a geometrical tilting effect, hydro- 
dynamic modulation effects being of second order. Provided the large-wave steep- 
ness Is small compared to the total surface roughness as measured by the total 
(diffraction-effective) mean square slope $ 2 , then the cross section variation 
will be proportional to the large-wave slope component In the plane of Incidence 
>rj »x as 
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where o is the angle of incidence and n° is the average cross section of the 
surface. In (?) the first term represents a linearized area tilt term while the 
second term represents the rigid rotation of the small-scale scattered power pat- 
tern by the large wave slopes. The average cross section is proportional to the 
probability density function of surface wave slopes evaluated at the specular 
condition for backscatter (namely, slope = tan o). Assuming a Gaussian isotropic 
distribution of slopes, the cross-section roll-off is given approximately by 

1 3o° = - 2 tan •) (3) 

7 ° 3 o a 2 

where again fl 2 is the total mean square slope. 

Assuming that the water wavelength 2ir/k is small compared to the azimuth 
besmwidth Ly, kLy » 1, it follows that the spectrum of m(x, 4 ) is proportional 
co the directional wave slope spectrum. If the gain pattern is assumed to be 
Gaussian, G(y) = expf-y 2 /2Ly ? - ) , one finds that in the limit of larqe kLy the dir- 
ectional modulation spectrum is given by 

p m( k » * ) = [cot r ’ - l_ tan _ e l 2 k 2 F(k, *) (4) 

Ly L G 7 * J 

where F is the polar-symmetric directional height spectrum, defined such that 
the height variance. 


<r 2> = A)kdk^ (S) 

We note that in (1) the azimuth coordinate y was treated as rectilinear; this Is 
permissible for directionally spread seas. The wave front curvature enters, along 
with the finite footprint size and antenna rotation, in determining the direc- 
tional resolution. In the satellite case the resolution Is typically 7° (200 m 
water wave), while in the aircraft case the resolution Is typically 20°. 

The ROWS data processing, described in detail in Jackson et al. (1984a), 
consists of first correcting for the wave front sphericity on a puTse by pulse 
basis on going from the signal delay time to the surface range coordinate x, and 
then integrating the pulse returns in surface-fixed range bins over a time cor- 
responding to 15° of antenna rotation (N * 42 pulses). The surface tracking Is 
accomplished usinq an input aircraft speed. The modulation m(x, 4) is then com- 
puted by normalizing by an estimate of the average power envelope <W^x, 4> ob- 
tained by averaging over several antenna rotations. Unity Is then subtracted and 
the data are rewindowed and the spectrum P m ( 1c , 4) computed by fast Fourier trans- 
form. The final estimates of P m are obtained by averaging the spectra over sev- 
eral antenna rotations, subtracting a computed residual fading spectrum, and cor- 
recting for the finite pulse response (20% spectral roll-off at 40 m wavelength). 

We note that the tilt model solution (4) has been shown by Jackson (1981) to 
correspond to the first term In a series expansion of the geometrical optics 
solution for P m , where the orderlnq parameter Is the large-wave steepness. Non- 
linear terms, both electromagnetic and hydrodynamic, were found to be small pro- 
vided that 1) f 2 Is sufficiently larqe (wind speeds greater than ca. 5 ms -1 ), 11) 
the Incidence anqle lies In the range of 8°-16°, and 111) the larqe-wave slope 
Is not too larqe (steepness < 0.10). The ROWS technique In general, and the 
tilt model result In particular, have been extensively validated by Jackson et 
al. (1984a). In the following, we will Illustrate the power of the ROWS technique 
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Our present understand nq of the Tasic physics of ocean waves is to a large 
extent embodied in numerical wave forecast/hindcast models, specifically, in the 
riqht hand side of the spectral energy transport equation. Considerable uncer- 
tainty exists in the parameterization of wave spectral growth (and decay) under 
the action of variable winds, particularly in the specification of directional 
distributions and the redistribution of wave nnerqv over direction in furninq 
wind fields. This is abundantly c.lnar f rom Hasselmann's (IQ34) model intercom- 
parison report. We have undertaken a study with V. Cardone "Iceanwpaf her , Inc.) 
comoarinq ROWS data taken in the Norwegian Sea during an intense storm to special 
hindcast runs with wel 1 -speci f ied wind fields (or as w *1 1 -spec i find as possible). 
The comparison is particularly interesting because the storm not only produced 
very high seas (ca. 10 m) , but created a complex, rapidly evolving wave field. 
The comparison thus provides a strong test of model performance. Initial model 
runs have been made with a fine-mos^ OHO km/"* hr) version of Cardone e* al.'s 
(1976) ODOP discrete spectral model (a variant of the U. S. Navy's opera t Tonal 
SOWM model) and with a coupled discrete model, the SAIL model (cf. Hasselmann, 
1984). The hindcasts used all available ship wind reports on record and aircraft 
winds obtained during a low-level flight leg. The modeled area was the N. E. 
Atlantic, and the spin up time (before the flight of interest) was two weeks. 

The NASA Ames' CV-990 flight track and ROWS data takes between ca. 0800 Z 
and 1000 Z on Nov. 3, 1978 are shown in Fiqure 2 along with the nearby 0DGP model 
grid points. The flight track constitutes a box pattern measuring ca. 150 km 
(N-S dimension) by 700 km (SW-NE direction). The aircraft altitude was ca. 10 km. 
Ten ROWS files were obtained in the flight box; these are indicated by the let- 
ters A-J in Fiqure P. At the northern end of the box, close to file A, is the 
weather station Tromsoflaket (TR0MS0) which provided 3-hourly wind reports and 
Waverider wave spectra observations. The synoptic situation at 1800 Z on Nov. 2 
at about the peak of the storm is shown in Figure 3. By flfqht time, the low 
had moved to the north of the box, and the winds had shifted to westerly over 
the northern end of the box. The high pressure ridge also moved northward. Its 
axis transectlnq the midsection of the box at the time of the fliqht. Figure 4 
Is a presentation of the ten files of ROWS data in the form of polar contour 
plots of directional height spectra laid out on a map of the Norwegian Sea. The 
spectra, contoured at equal intervals, are scaled to the peak values; absolute 
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energy levels are indicated by the estimated significant wave heights given 
numerically in meters beside each spectrum. The spectra shown are one-sided. 
The ISO 0 ambiguity in the ROWS spectra was resolved and the true direction of 
wave travel determined on the basis of the synoptic situation. The synoptic 

view of the wave field in Figure 4 shows two major wave trains traveling at 
nearly right angles to each other. The NE travelling component appears to enter 
the box at the southern end more or less as swell, decay slightly, and then qrow 

rapidly in the stronqer winds in the northern end of the box, where the winds 

are westerly at ca. 20 ms' 1 . The NF direction of this wave train is maintained 
over the entire 700 km distance of the flight legs. Only at the northernmost 
location does one see energy appearing in the local wind direction (this is 

evident in the comma shape of file A). The SE traveling component is evidently 

a pulse of wave energy generated in the western sector of the rapidly moving 
cyclone at about the 1800 Z map time (Figure 3). 



I 

l. ROWS Directional Heiqht Spectra in Norwegian Sea Storm, Tic car' , i 

are 0.05 Hz. Spectra are scaled to peak values. The numbers 1 r. - ! 

dicate wave heights in meters. 
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An indication of the quality of the ROWS data is given by Figure 5, which 
compares the nondi rectional spectrum of file A/subfile a with the TROMSO Wave- 
rider spectrum. The agreement is seen to be excellent. 


The results of the first ODGP model run showed basically good agreement 
with the ROWS and TROMSO significant wave heights and nondi rectional spectra. 
However, comparison of directonal spectra showed immediately that the model had 
failed to produce the SE travelling component in the strength observed. This was 
seen to be due a poor kine'atic analysis for the data-sparse western sector of 
the cyclone. Figure 6 compares the ROWS spectrum for file b with the ODGP hind- 
cast spectrum for grid point 254 (0900 Z) in FNOC/SOWM/ODGP variance format. 
Only a trace of the SE travelling component is seen at higher frequency. We see 
further in Figure 6 that the hindcaot is also severely underestimating the 
strength of the NE travelling component, which is only weakly represented in the 
hindcast as an islet of energy near 0.065 Hz. The hindcast is evidently putting 
the energy of the NE component (to 45°) into the local wind direction (to 90°). 
While the model errors here are ob/ious, it is important to note that they are 
n^ten masked in the nondirectional spectrum: Thus, the nondirectional spectra 
corresponding to Figure 6 shown in Figure 7 are quite close and give no clue as 
to the seriousness of the actual hindcast error. That the model should perform 
better for the integrated properties should not be surprising since the greatest 
uncertainty in wave models today lies in the modelling of the directional re- 
sponse to varying winds (the directional relaxation problem). 


A second ODGP model run with a new kinematic analysis of the western sector 
of the cyclone essentially correctly reproduced the observed SE directed compo- 
nents. However, the model still tended to place more energy in the local wind 
direction than actually observed. A first run of the SAIL model exaggerated 
this tendency in the ODGP. Additional runs with the SAIL model will be made 
with slower directional relaxation rates; also the the nondirectional growth 
rates will be tuned to maximize the agreement with the ROWS (and TROMSO) wave 
observations. In sumina r y. this hindcast comparison shows i) how the ROWS obser- 
vations pointed immediately to wind field specification errors, and ii) that the 
directional relaxation rates, especially of the newer-generation coupled models 
(cf., Hasselmann, 1 984 ) are too fast. This latter observation is supported by 
the results of the next case study. 



Figure 5. Comparison of ROWS Non- 
directional Spectrum 
with TROMSO Waverider. 
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FNOC FORMAT VARIANCE SPECTRA HINDCAST VS. RADAR SPECTRUM 



Fiqure 6 . Figure 7 - 

3. FETCH LIMITED WAVE GROWTH ORSERVED DURING MASEX 

The OSFC rvdar was installed on the NASA Wallops' P-3 aircraft in 198? with 
an improved am.enna system installation (which eliminated diffraction problems 
experienced in the Fall '78 Mission data). The main reason for the switch of 
platforms for the ROWS was to be able to fly in concert with E. Walsh's Surface 
Contour Radar (SOI), a 36 GHz, direct topographic mapping radar, in order to ob- 
tain high resolution directional spectra for intercomparison/va i idation purposes. 
The ROWS participated in three joint flights with the SCR during the MASEX (Meso- 
scale Air-Sea Exchange) experiment in January 1983. Besides systems intercompar- 
ison, the flight objective was to obtain data on the evolution of the directional 
spectrum with fetch Coring strong cold air outbreaks off the east coast of the 
U. S. Figure 8 shows the ROWS flight track (7 km altitude) for the 1/16 MASEX 
flight, the tape/file numbers being indicated by numbers and letters. The wind 
was approximately 12 ms -1 blowing offshore normal to the coast and to within 
a few degrees of the flight track. Figure 9 is a selection of the ROWS direction- 
al height spectra from the 1/16 flight. Figure 10 shows the entire sequence of 
-7 ROWS-inferred nondirectional spectra for the flight leg extending out to c a. 300 

km. The ROWS directional spectra of Figure 9 show two nearly equally energetic 
t wave components, one travelling downwind, the other travelling in directions 

nearly opposite to the 1 ine-of-sight directions to the mouth of the Delaware Bay. 
The angles of the 'Delaware Bay' component are seen to be steeper than the llne- 
of-siqht angles to the mouth of the bay; this is apparently due to refraction 

j near the mouth of the Bay. The existence of this strong off-wind component was 

-J first discovered by Walsh <}t. al. (1982). It Is evidently due to fact that the 

.j waves outside the mouth of the~Fav have a 'leq up' over the waves further south 

.3 alonq the coast. These waves may then enter a more rapid growth stage (Miles' 

5 qrowth) earlier than che downwind waves along the track. There Is also the pos- 

^ Iblllty that these waves are preferentially qrown (In file 4-c these ..aves are 

-f 
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Figure 8. Flight Track and ROWS 
Data Takes for MASEX 
1/16/83 Flight. Wind 
normal to coast anj 
parallel to flight 
track. 


seen to exceed the downwind component) because they are very close to the 
Phillips' resonance curve (shown as the horizontal line in Fiqure 9). Most im- 
importar.t, we observe in these spectra that over the entire fetch the off-wind 
component does not turn into the wind direction: It is not ’absorbed' by the 
'local wind sea' (downwind component) until its mean direction lies well within 
the directional spread of the downwind component. This observation runs counter 
to current theory (Hasselmann, 1984) which would ascribe a strong directional 
coupling to the two components. On the contrary, the two components appear to 
be uncoupled, or at the most only weakly coupled. This behavior is consistent 
with what we have observed in the Norwegian Sea hindcast study. Preexisting 
components are not easily turned into the wind direction, nor do they tend to be 
absorbed by the local wind sea. This conclusion was reached also by Holthuijsen 
(1983) who found strong off-wind components in fetch-limited spectra in the North 
Sea, also due to coastline irregularities. The ROWS (and SCR) spectra from the 
two other MASEX flights similarly show strong off-wind components that point 
clearly to major embayments in the Middle Atlantic Bight. Similarly, also, 
these components appear to be decoupled from the downwind component. Figure 11 
is a ROWS height spectrum from a flight down-fetch of Long Island which shows 
dominant wave energy not it, the wind direction (to 160°), but at an angle point- 
ing to the center of the Block Island/Rhode Island Sound complex. 

The ROWS data of Figure 11 can be compared to nearly colocated SCR data in 
Figure 1?. The overall agreement is seen to be excellent. Slight discrepancies 
are most likely due colocation discrepancies, both spatial and temporal. It is 
worth noting here that apart from stereo photoqraphy, only the SfR could have 
provided this kind of high resolution intercomparison data. 
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Figure 10. ROWS Nondi recti onal Spectra from 
1/16 MASEX Flight for Unit Sen- 
sitivity Coefficient. 



Fiqure 9. Examples of ROWS Directional Height 
Spectra f rom 1/i ^3 MASEX Flight. 

From top to hot' :he files shown 
are (cf. Figure 8/ 5B, 4C, and 3C. 
The aircraft heading, equals the 
(anti) wind direction, is to the top 
of the page. The arrows indicate 
the 1 ine-of-sight anqles to the 
mouth of the Delaware Bay. The 
horizontal lines are the Phillips' 
resonance curve. The spectra are 
contoured with six equally spaced 
contour Intervals scaled to the peak 
values. Frequency rings are 
0.10 Hz and 0.20 Hz. 
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Figure 11. ROWS Height Spectrum from 
1/20 MASEX Flight. 81 C is 
Block Island Component. 


Figure 12. SCR Height Spectrum for 
Comparison to Figure 11. 
Colocation is not exact. 


A thorough analysis of the MASEX data set should go a long way to resolving 
some fundamental issues in the physics of wind-generated ocean waves, at least 
in the macroscopic, or practical modelling sense. 

4. 'THE SON OF SWOP 1 — AN OBSERVATION OF A FULLY DEVELOPED SEA 

A detailed analysis of a Fall '78 Mission ROWS observation of a fully de- 
veloped sea in the N. E. Pacific and intercomparison with a pitch-roll buoy 
is given by Jackson (1984b). The ROWS-observed directional spectrum of the 3.3 
m significant wave height sea is given in Figure 13, and the nondirectional com- 
parison with the buoy is given in Figure 14. The directional comparison (not 
shown here) gave excellent agreement in terms of the mean wave directions and 
directional spreads as functions of frequency. 

The observed spectrum in Figure 13 bears such a remarkable resemblance to 
the classical SWOP (Stereo Wave Observation Project) spectrum (Pierson, 1960) 
that we have dubbed this ROWS spectrum 'The Son of SWOP'. The two spectra have 
nearly Identical half power widths (80° at peak to 120° at high frequency), and 
both spectra exhibit nearly identical bimodal structures which accord with the 
Phillips' resonance condition. The Phillips' resonance angles y , determined 
according to the resonance condition, 

U c cos y = c(f) = q/2*f (6) 

where U c is the windspeed, or convection velocity, and c is the water wave phase 
speed, are Indicated by the arrows in Figure 15, where U c * 14.2 ms" 1 is chosen to 
<j1v3 y » 0 at f ■ 0.11 Hz. The agreement with the observed modal angles Is seen 
to be very good. This observation Is seen to confirm Phillips' original conten- 
tion (Phillips, 1958) that the SWOP spectrum bimodality was a real manifestation 
of the resonance mechanism and not a statistical fluke. 
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Fiqure 13. 'Son of SWOP' Height Spectrum. Figure 14. 
Contouring as in Figure 9. 


Nondi recti onal Heiqht 
Spectrum Corresponding 
to Fiqure 13. Solid line 
is radar, circles buoy. 
Dots are radar confidence 
intervals. 




Figure 15. 'Son of SWOP' Directional Clits 
(Normalized to Peak) 
with Philips' Resonance angles 
Indicated by arrows. 
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5. CONCLUSION 


This paper has emphasized the use of the ROWS in case studies of wave physi- 
cal processes. The results, while preliminary, give an idea of the enormous 
potential of this remote sensing technique. Further study of the Norwegian Sea 
case and the MASEX data shoud permit more concrete, quantitative statements 
regarding wave growth and directional relaxation. The 'Son of SWOP' observation 
fairly conclusively demonstrates that the Philips' resonance mechanism is indeed 
operative in the the fully developed sea state and is a major factor in determin- 
ing the directional structure. 

We have not dwelt on the details of the ROWS technique here, its accuracy, 
or its ultimate space application. These aspects have been extensively delt 
with in Jackson (1981) and Jackson et jsK (1984a). However, the reader may 
appreciate the remarkable accuracy of the indirect ROWS measurements by again 
examininq the Waverider buoy, SCR, and pitch-roll buoy comparisons given respect- 
ively in Fiqures 5, 12, and 14: The technique appears to be very accurate, at 
least for sea states above 1 m, and providing the windspeed is greater than ca. 
5 ms -1 . Similar accuracy is possible with a spacecraft system employing present 
generation short pulse radar altimeters. 
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